Nematicity from mixed S± + d x 2_ 2 states in iron-based superconductors 
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We demonstrate that in iron-based superconductors (SC), the extended S± SC state coexists with 
the d x 2_ y 2 state under generic conditions. The mixed S± + d x 2_ y 2 SC is a natural nematic state in 
which the tetragonal symmetry C4 is broken to C2 explaining puzzling findings of nematic SC in FeSe 
films [l|]. Moreover, we report the possibility of a first order transition at low-T from the nematic 
S± + d x 2_ y 2 state to the pure d x 2_ y 2 state induced by the Zeeman magnetic field proposing an 
original experimental strategy for identifying our mixed nematic state in FeSe films. Extrapolating 
our findings, we argue that nematicity in non superconducting states of underdoped and undoped 
pnictides may reflect mixed S± + d x 2_ y 2 Density Wave states. 
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Iron-based superconductors (SC) support high criti- 
cal temperatures and unconventional superconducting 
states like cuprates. However, while in high-T c hole 
doped cuprates the gap is certainly d-wave, in iron- 
based materials the situation is much more complex 
and still under investigation. In a series of materials 
like Bai_ :E K :E Fe2As2, BaFe2- x Co :c As2, FeTei-^Se^ and 
K x Fe2-^Se2 there is experimental evidence for nodeless 
superconductivity Q. Apparently, this is not the usual 
isotropic s-wave gap but it is instead the extended s- 
wave S ex t = S± gap changing sign between the electron 
and hole Fermi surfaces [3j. On the other hand, in ma- 
terials like KFe 2 As 2 , LiFeP, LaOFeP, BaFe(Asi_ ;r P :c )2 
and BaFe2- ;E Ru :c As2 experiments as diverse as NMR, 
STM, thermal conductivity and penetration depth mea- 
surements all point to the presence of gap nodes on 
the FS mill. On the theoretical side, unconventional 
either nodeless or nodal gap structures are usually re- 
garded as evidence of non-phononic mechanisms 0, HI, • 
Nonetheless, we have reported recently [7] that electron- 
phonon interaction dominated by small-q processes pro- 
duces nodeless S± as well as nodal SC states depending 
on the doping, including a phonon- driven triplet p-wave 
state 0] possibly observed recently in LiFeAs 0|. 

Among the iron-based SC, FeSe is the simplest com- 
pound representing a prototype for this class of materials 
on which ideas may be tested. Recent scanning tunnel- 
ing microscopy and spectroscopy measurements on high 
quality FeSe films by Song et al. [1] reveal an aston- 
ishing feature. Apparently, the tetragonal symmetry C4 
is broken to C2 pointing to a nematic SC state, that 
can certainly not be attributed to the tiny orthorhom- 
bic distortion of the FeSe lattice. Nematicity in iron- 
pnictides is already a highly debated issue following re- 
orts for a possible electronic nematic phase transition 
that often coincides and possibly drives [10-12] the 
orthorhombic distortion which accompanies the antifer- 
romagnetic transition in undoped and underdoped iron 
pnictides. Because of the coincidence of nematicity and 
antiferromagnetism, it has been suggested that the mag- 
netism itself drives an electronic nematic phase transition 
[l^4l(|. However, in the case of FeSe there are no such 



antiferromagnetic phases involved. As a matter of fact, 
there have been also proposals for antiferromagnetism- 
independent-nematicity such as ferro-orbital nematic or- 
dering [171 ]. 

In the present Letter we introduce a novel approach to 
the phenomenon of nematicity in pnictides. We demon- 
strate that S± SC coexists under generic conditions with 
d x 2_ y 2 SC and the mixed S = j-+d ;r 2_ 2y 2 SC state is a promi- 
nent nematic state capable of explaining the reports of 
nematic SC in FeSe films [1]. Moreover, we point out 
that at low temperatures, a Zeeman field can induce a 
first order transition from the nematic $±+d x 2_ y 2 state 
to the pure d x 2_ y 2 state. We therefore propose an ex- 
perimental strategy for identifying the mixed S±+d x 2_ y 2 
character of the nematic SC state by repeating the ex- 
periments of Song et al [1] in the presence of in-plane 
magnetic fields in which case the Zeeman effect will dom- 
inate if the FeSe films are sufficiently thin [l8j. Quite 
remarkably, we were able to produce self-consistently a 
Zeeman Field - Temperature phase diagram associated 
with the mixed nematic phase that exhibits three dis- 
tinct SC regions and a tetracritical point reminding the 
one observed in UPt3 [19]. Finally, we show that singlet 
mixed SC states that violate time-reversal symmetry (T) 
are also accessible, exhibiting Zeeman field- induced tran- 
sitions within the SC phase. 

For our generic discussion, we can model qualita- 
tively iron-based SC like FeSe with a minimal and suf- 
ficient two-band model exhibiting a hole pocket around 
the r(0,0)-point and an electron pocket around M(7r, tt)~ 
point: s e (k) = ti(cos & x + cos k y ) — £2 cos k x cos k y -\-C — \i 
and Eh{k) = t\ (cos k x + cos k y ) + £2 cos k x cos k y — C — fi 
where \i is the chemical potential and we set t\ = 1, 
t2 = 0.5, C = 2. The present dispersions capture the 
necessary ingredients allowing for the S± state to emerge, 
namely well separated electron and hole Fermi surface 
(FS) sheets. Moreover, we have verified that the results 
that we report here are independent of any further band 
structure details. In fact, we keep our analysis as generic 
as possible by adopting a separable potentials approach 
for the effective interactions, allowing for a broad discus- 
sion concerning a large number of iron-based SC. 
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For this two band system we can write 



H = [ee(fc)4 >ff c fciff + e h {k)dlj k ^] - | £ V(fe, ft') (4 



k,a 



fc,fc' 



(i) 



where } a and dj^ a are annihilation (creation) operators 
for the electron e e (k) and hole band Sh(k) bands respec- 
tively of spin projection a =t,i and N the number of 
lattice points. Notice that for the specific choice of the 
interaction [20|, where intraband and interband interac- 
tion strengths are equal, electron and hole bands share 
the same superconducting order parameter A(fe), which 
is defined in the following manner 

A ( fe ) = -^Y. v ^ k '){ c -y,i c k'A + d -y,i d y,t) 

k' 

A *( fe ) = -^E 7 ( fe > fc ')(4A4+4A^(2) 

k' 

The order parameter A(fe) = ^ n A n f n (k) consists of ir- 
reducible representations (IR) f n {k) of the point group 
D4/i, that includes C4 as a subgroup, and is appropriate 
for describing the normal (tetragonal) phase of strongly 
two-dimensional pnictide compounds. Here we shall re- 
strict our analysis to the following IRs: fs(k) = 1 
(Aip), fs±(k) = cosk x + cos k y (A lg ) and fd x2 _ y 2( k ) = 
cosk x — cosky (Big). All of them are even under in- 
version as it is required for singlet intraband supercon- 
ductivity and connect up to nearest neighbors. Note 
that by Fourier transforming the corresponding effec- 
tive interaction field V(k,k') = ^Z n V n f n (k)f n (k f ) with 
n = 5, S±, d x 2_ y 2 one can see that in the context of 
real space extended Hubbard models this would corre- 
spond to on-site interactions of the form Vs = U 5^ and 
nearest neighbor interactions of the form Vs±, Vd x2 _ y2 = 

Vi 5 j(^j^ + (±i ? o) +^j,i+(o,±i))j where i and j the real-space 
lattice points indices. Nevertheless extended Hubbard 
models are not the unique option since also a small-q 
phonon-mediated pairing potential leads to a plethora 
of non s-wave IRs such as S± or d x 2_ y 2 0, [22(- With 
this reasoning we conclude that by studying here the 
SC phase competition via separable potentials we report 
generic results that are universally relevant, and indepen- 
dent of the exact microscopic mechanism of SC. 

Within the aforementioned subspace of IRs, there are 
only two minimal schemes to achieve a nematic state. 
These are the mixed states S + d x 2_ y 2 and S± + d x 2_ y 2. 
Of course, the cases S + S± + d x 2_ y 2 1 iS + S± + d x 2_ y 2, 
S + iS± + d x 2_ y 2 are also possible but not minimal. In 
all these symmetry breaking patterns the subgroup C4 
reduces to C2. Notice that for a minimal nematic phase, 
the two IRs involved must lock in the same phase. If the 
two phases lock in phases with it/ 2 difference then the 



mixed state leads to broken T but unbroken C4. These 
states are also important and for completeness we shall 
also discuss features of their phase diagram and their 
phenomenology. 

We consider first the minimal configurations in which 
nematicity emerges and T is preserved. In these cases 
the order parameter A(fc) is the sum of two order pa- 
rameters Ai ? 2(fe) corresponding to the two coexisting 
IRs. For simplicity we shall consider that Ai^fe) are 
real. At this point, we introduce the spinor Sfr\, = 

( C L t ^fe t c — k I ^—k\) an d employ the usual Pauli ma- 
trices f , p. The mean field Hamiltonian can be written 
in the following compact form 

k v x ' 

-MT 3 po - #T po + Ai(fc)fip + A 2 (fc)fipo|^fc (3) 

where we also incorporated the effect of a Zeeman field 
B. With usual techniques we calculate Green's functions 
that exhibit four quasiparticle branches and coupled self- 
consistent equations that provide the two gaps Ai^fc). 

To illustrate the fact that achieving mixed states is 
not trivial, we start with the competition between the 
isotropic S IR and the d x 2_ y 2 IR. As expected, s and 
d-wave SC phases are highly competitive and coexistence 




Temperature 

FIG. 1: (color online) Typical temperature induced first or- 
der transition from S (red) to d x 2_ y 2 SC (green) for various 
chemical potentials fi obtained when Vs < Vd x2 _ y2 < l-5Vs. 
Free energy calculations not reported here verify the reality 
of this transition. 



3 



















































































































































































S 


-ext + d 2 2 
























y 











































































2 3 
V, 




I 





S-ext 



FIG. 2: (color online) a. Map of the phases obtained varying 
Vs± = Vs-ext and Vd x2 2 f° r zero temperature in the hole 
doped (fi = —0.4) regime (S ex t = S± (red), S±-\- d x i_ y i (yel- 
low), d x 2_ y 2 (green)). Notice the extended region over which 
S± and d x 2_ y 2 coexist, b. An example of self-consistently 
obtained total SC gap amplitude in the nematic S± + d x 2_ y 2 
state plotted on the electron and hole Fermi surface sheets for 
a hole doped case \i = —0.4. Note that the tetragonal symme- 
try C4 has been reduced to C2. c. Corresponding free energy 
results as a function of the d x 2_ y 2 and S± gaps (Ad x2 _ 2 
and As ± respectively) exhibiting four degenerate total min- 
ima for which both order parameters are finite confirming the 
S± + d x 2_ y 2 state. 



cannot be achieved at any value of the respective poten- 
tials. Remarkably, when Vs < Vd x2 _ y2 < 1.5 Kg we ob- 
serve a first order transition from d-wave to s-wave SC as 
we lower the temperature (Fig. 1). We confirmed the re- 
ality of this transition, as we did for all results reported 
here, by verifying that it minimizes the corresponding 
free energy. 

While the isotropic s-wave SC state cannot coexist 



FIG. 3: (color online) Self-consistently obtained nodal ne- 
matic S± + d x 2_ y 2 state (plotted over the whole Brillouin 
zone) using a fully momentum dependent small-q electron- 
phonon pairing potential and a realistic four band model for 
hole doped iron pnictides described in Ref. [3. 



with the d-wave state, the extended s-wave S± state that 
is widely considered relevant for iron-based superconduc- 
tors coexists with d x i_ y i over a wide range of the effec- 
tive potentials (Fig. 2a). We report in Fig. 2b a typical 
solution in the mixed S± + d x i_ y i state plotted solely 
on the Fermi surface for a hole doped system \i = —0.4. 
The emergence of nematicity is directly evident from Fig. 
2b. In Fig. 2c we present free energy calculations, ex- 
hibiting four degenerate (because of symmetry) minima 
corresponding to the mixed nematic SC state. We in- 
sist that all effective potentials and dispersions used in 
our self-consistent calculations 'preserve tetragonal sym- 
metry. Only because S± and d x 2_ y 2 coexist, tetragonal 
symmetry is broken and nematicity emerges. The essen- 
tial ingredients leading to the nematic S± +d x 2_ y 2 state 
is on one hand the well separated electron and hole pock- 
ets that favor the stabilization of the S± IR and on the 
other, some weak tendency towards the formation of the 
d-wave that is further assisted by the presence of S±. 
The detailed characteristics of the Fermi surface topol- 
ogy are not crucial for the formation of the nematic state 
but mainly determine the exact balance of the S± and 
d x 2_ y 2 OPs, that controls the nodal or nodeless type of 
the quasiparticle excitation spectrum. 

To demonstrate that our qualitative findings are not 
peculiar to details of the considered dispersion or to the 
separable potentials character of the above analysis, we 
also report in Fig. 3 a self-consistently obtained nodal 
nematic S± + d x 2_ y 2 state using a fully momentum de- 
pendent small-q phonon-mediated interaction V(k, k') = 
V£ b ~ V P h(q 2 c +\k- fef)" 1 (V Cb * 0.1V ph and q c = tt/6) 
and an accurate hole doped (/i = —0.6) four band model 
for high-T c iron pnictides described in [7[ . Naturally, also 
in this case both the dispersion and the interaction used 
in our calculations preserve the tetragonal symmetry, and 
only the resulting self-consistent solution depicted in Fig. 
3 exhibits nematicity. Note that small-q phonon driven 
unconventional SC has also been considered in the past 
for high-T c cuprates [22|, [23[ , heavy fermion [24| organic 




FIG. 4: (color online) Zeeman field - Temperature phase di- 
agram resulting from our self-consistent calculations for the 
case shown in Fig. 2b. It exhibits three distinct supercon- 
ducting regions and a tetracritical point as it is schematically 
depicted in the inset. A similar diagram has been experimen- 
tally observed in UPt 3 0. 

(25| and cobaltite [26| SC and is known to produce a 
loss of rigidity of the gap function in momentum space 
called momentum decoupling 22} thus allowing for gap 
symmetry transitions. Note also that depending on the 
relative magnitude of the two order parameters, which in 
turn depends on the effective potentials, the mixed ne- 
matic state S± + d x i_ y i can be either nodeless as in the 
example of Fig. 2b or nodal as in Fig. 3. The small- 
q results are just a particular case confirming that the 
findings of the separable potentials analysis is generic. 
The nematic S± + d x 2_ y 2 SC state is indeed a model in- 
dependent phenomenon likely to be behind the puzzling 
reports of nematicity in FeSe films 

Moreover, the qualitative behavior of the mixed S± + 
d x 2_ y 2 state in the presence of a Zeeman field may al- 
low its firm identification from the experiments. At low 
temperatures, we obtain a remarkable first- order field- 
induced transition from the nematic S± -\-d x 2_ y 2 state to 
the pure d x 2_ y 2 state (Fig. 4). The transition exists only 
at sufficiently low temperatures. We therefore propose an 
experimental approach for identifying the S± + d x 2_ y 2 
state by applying in-plane magnetic fields to the FeSe 
films instead of perpendicular fields as in Exploring 
the higher temperatures in the presence of the field al- 
lows to construct with self- consistent calculations Field 
- Temperature phase diagrams. Quite remarkably, for 
the zero field state shown in Fig. 2b, we obtain a Field- 
Temperature phase diagram exhibiting three distinct SC 
regions and a tetracritical point (Fig. 4) in analogy to 
the well known phase diagram of UPt3 obtained there as 
well in the presence of in-plane fields [l9|. In our case 
there is a simple understanding of this diagram as fol- 




FIG. 5: (color online) Field and temperature behavior of 
the mixed S + id x 2_ y 2 state obtained for Vd x2 2 sufficiently 
stronger than Vs- Such a mixed state breaks time reversal in- 
variance but preserves the tetragonal symmetry. Here as well 
there is at low-T a first order Zeeman field induced transition 
from the mixed to the pure d x 2_ y 2 state. The inset depicts 
schematically the corresponding phase diagram. 



lows: The S± state is stronger in this example owing the 
higher T c at zero field. On the other hand, in the pres- 
ence of a large Zeeman field, the d x 2_ y 2 state with nodes 
on the FS is energetically more favorable than a nodeless 
SC state [25], like S±, which exhibits a higher critical 
field at zero temperature. Therefore, the reason for such 
a complicated phase diagram lies in the extraordinary 
fact that at zero field, a fully gapped state like S± allows 
at lower temperatures its coexistence with an emergent 
nodal d-wave state. 

Mixed states may also lead to T breaking despite the 
singlet character of the condensates. This naturally hap- 
pens if the two coexisting order parameters lock in a tt/2 
phase difference. As we have already mentioned, when 
two different IRs coexist with this type of phase locking 
C4 symmetry is preserved. The mean field Hamiltonian 
describing such a situation has the following form 

k v x ' 

-MT 3 po - Sf po + Ai(fc)fip - A 2 (/c)f 2 po|^fc (4) 

leading to a different system of coupled self-consistent 
gap equations, compared to the previously examined 
case. We briefly report some results in Fig. 5 that will 
help the interpretation of eventually observed in-plane 
field-induced first order transitions in the experiments 
that we suggest in the previous paragraph. While the 
state S + d x 2_ y 2 is not accessible, breaking of T due 
to 7r/2 phase locking, allows for sufficiently large d-wave 
potentials (Vd 2 _ 2 > l-2Vs) self-consistent solutions of 
the S + id x 2_ y 2 (or equivalently iS + d x 2_ y 2) type. Re- 
markably, here as well, we obtain at zero temperature 
a first order Zeeman field-induced transition from the 
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mixed S-\-id x 2_ y 2 (iS + d x i_ y i) state to the pure id x i_ y i 
(d x 2_ y 2) state (Fig 5) that exhibits a qualitatively sim- 
ilar behavior to the field-induced transition from the 
S± + d x 2_ y 2 state to the pure d x 2_ y 2 discussed in the 
previous paragraph. Nevertheless, the mixed S + id x 2_ y 2 
state preserves the tetragonal symmetry and therefore a 
first order Zeeman field-induced transition within the SC 
phase does not necessarily imply the presence of a ne- 
matic SC state. The Zeeman field - temperature phase 
diagram depicted in Fig. 5 exhibits now only two distinct 
SC regions because the nodal d x 2_ y 2 state with the higher 
critical field is now stronger at zero field than the S± 
state having the higher critical temperature as well. On 
the other hand, for the competition of d x 2_ y 2 (id x 2_ y 2) 
with iS± (S±) we observe no qualitative difference in the 
phase diagrams compared to the one obtained previously 
when the two order parameters lock in the same phase. 
In fact, it has been suggested that enhanced pnictogen 
height may favor this type of T violating mixed states 

fir 

Finally, from the point of view of BCS theory, un- 
conventional particle-hole condensates like spin singlet 
or triplet Density Waves (DW) behave similarly to SC 
condensates when the two bands are perfectly nested. 
Qualitatively, our present findings could be extrapolated 
to the DW condensates suggesting that the nematicity 
that accompanies the antiferromagnetic transition in un- 
doped and underdoped iron-pnictides may well indicate 
the presence of a nematic mixed S± + d x 2_ y 2 spin DW 
state. Note that in particle-hole asymmetric systems a 
spin DW and charge DW exhibiting the same momentum 
structure coexist [27j , and a nematic charge DW may in- 
deed drive the orthorhombic distortion in the form of a 
Peierls instability [10]. If it is proven that unconventional 



SC in iron based SC emerges in the proximity of un- 
conventional density wave phases of the same symmetry, 
then a fundamental analogy to high-T c cuprates emerges 
where there are reports of a d-wave DW (orbital antifer- 
romagnet) associated with the pseudogap (28j |. Uncon- 
ventional density wave states may host exotic phenom- 
ena like emergent chiral density wave stat es [ 29j induced 
by the orbital effects of a magnetic field [30| producing 
extraordinary phenomena like the anomalous Nernst sig- 
nal (3l| observed in underdoped cuprates [32]. Anoma- 
lies in the Nerst signal are apparently present in under- 
doped pnictides as well at temperatures where nematicity 
emerges [33|. Dedicated work is needed in order to sub- 
stantiate the conjecture of a mixed nematic density waves 
state in pnictides. 

In conclusion, we claim that nematicity in iron-based 
superconductors indicates the presence of mixed conden- 
sates in which d x 2_ y 2 and S± order parameters coex- 
ist. We suggest the experimental search for an in-plane 
field-induced melting of nematicity in films of FeSe by 
an abrupt first order transition at low temperatures. 
Mixed states that break time reversal invariance are also 
shown to be accessible and should be taken into con- 
sideration in the analysis of experiments. If our find- 
ings for nematic mixed SC states are proven to be rel- 
evant in FeSe films, then it is probable that nematicity 
in non-superconducting phases of underdoped pnictides 
originates from analogous mixed nematic density wave 
condensates. 
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